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J  The  ^e*eent-ly-proposeci?algorithm  of  Karmarkar  -flj  for  solving  linear 
programs  ueee— the -following  step:  suppose  a  *  (a^,«..,a  )  >*)>  is  an 
interior  point  to  the  feasible  region  of  the— following  linear  program  Pr^ 

A 


minimize 
subject  to 


t 

C  X 


Ax  **  0 


x  >  0 


A  projective  transformation  maps  the  point  a  to  the  center  of  the 

* 

simplex  a_  and  leads  to -the~-fol  lowing  nonlinear  programs 

”  A  V 


TP: 


minimize  f(x')  *  c  ***- 

e  Dx' 

subject  to  ADx'  =  0 

eSc'  *  1 
x'  >  0 


where  D  is  a  diagonal  matrix  with  d^  *  a^.  TP  is  approximated  by  the 
following  program: 


AP: 


minimize  cCDx' 
subject  to  ADx'  ■  0 
eCx'  -  1 
x'  >  0 
lx'-aQl  <  or 
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Note  that  6  occurs  at  a  »  0.2451.  Further,  for  a  *  1/4  we  have  6  > 
max 

1/8  so  that  f(3g)  -  f(ag-arCp)  >  1/8.  However  one  might  note  that  f(x') 
cCDx' 

*  -  is  quasiconcave  over  the  positive  orthant  and  the  feasible  region 

eCDx' 

of  (TP)  is  convex  and  compact.  Hence,  the  optimal  solution  to  (TP)  is  a 
boundary  point.  Consequently  it  seems  logical  that  we  should  solve  AP 
for  the  maximum  possible  value  of  a  gaining  a  greater  reduction  of  the 
objective  function. 


Note  that  even  though  6  gives  a  good  approximating  bound  for  the 
decrease  in  the  objective  function  for  a  <  1/4,  the  bound  for  a  ■  1/4  is 
valid  for  all  a  >  1/4.  Hence,  the  order  of  complexity  is  not  adversely 
affected. 


A  second  logical  approach  is  to  use  the  maximum  value  of  a  but  to 


proceed  along  the  negative  gradient  direction  of  the  objective  function 
in  TP.  It  seems  probable  that  this  will  yield  faster  convergence  even 
though  it  is  not  guaranteed  to  do  so. 

EXAMPLE 

To  illustrate  the  points  made  above  we  solved  a  simple  problem 
proposed  by  Charnes  et  al.  [2]  using  each  of  the  above  approaches  and 
with  different  values  of  a. 

Here,  the  original  problem  is 


P:  minimize  x ^  s.t.  x^  +  X£  +  “  1  x^  >  0  i  ■  1,2,3  and  under 

a  projective  transformation  T  which  maps  an  interior  point  a  - 

cl 

(a^.a^ja^)  of  the  feasible  region  of  P  to  the  center  of  the  simplex  S 

-  {(y1,y2»y3):  yL  +  y2  +  y3  "  1*  y±  >  0  for  1  *  1>2»3>  we  set 


TP: 


minimize 


a2y2 


alyl+a2y2+a3y3 


s.t.  yeS 


Alternative  1:  Proceed  along  the  projected  negative  gradient  -c 
Following  Karmarkar's  notation,  we  have  B  »  [1,1,1]  so  that 


of  AP: 

P  * 


ri-Bt(BBt)  ^B]  DC.  Letting  c  ■  c  /lc  I,  and  noting  that  the  current 

P  P  P 

point  is  b  -  T  (a)  -  (1/3 ,1/3 ,1/3) C,  we  have  the  new  point 
d 


y  -  b  -  arc 

P 


(—  +  —  1/3  -  — 
'3  6’  '  3 


t 


where  0  <  a  <  1  and  r  -  l//n(n-l)  ■  1//F.  In  the  original  space  this 


point  is 


The  corresponding  point  in  the  original  space  is  then 


Dy 


(—7+  3ara,a./lld  It)/A 
3  1  2  g 

a2 

(— |  +  Sara^a^+a^)/  Id g>/A 


(— r  +  3ara„a_/»d  ll)/A 
3  2  3  g 


-  t  1  2  2 

where  A  =  e  Dy  =  -j  +  Oara^/  Hd^ll )  ^-^(a^a^)  +  a^]»  The  computational 

results  for  the  two  alternatives  are  shown  in  Tables  1  and  2  for  a  *  0.1, 

0.25,  0.50,  0.80,  and  0.90.  The  entries  in  the  tables  are  the  objective 

function  values.  The  starting  point  in  each  case  is  (0.10,  0.30,  O-SO)1". 

This  simple  example  strongly  supports  the  use  of  a  value  of  the 

A 

order  of  0.90.  The  results  using  the  direction  -d  is  comparable  with 

& 

A 

that  obtained  using  -cp*  Clearly  further  computational  testing  is  needed 
before  any  firm  conclusions  can  be  drawn.  The  order  of  complexity  under 
alternative  2  also  needs  to  be  investigated. 
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